Projected wavefunction study of Spin-1/2 Heisenberg model on the Kagome lattice 
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We perform a Gutzwiller projected wavefunction study for the spin-1/2 Heisenberg model on the 
Kagome lattice to compare energies of several spin-liquid states. The result indicates that a U(l)- 
Dirac spin-liquid state has the lowest energy. Furthermore, even without variational parameters, 
the energy turns out to be very close to that found by exact diagonalization. We show that such a 
U(l)-Dirac state represents a quantum phase whose low energy physics is governed by four flavors 
of two-component Dirac fermions coupled to a U(l) gauge field. These results are discussed in the 
context of recent experiments on ZnCu3(OH)6Cl2. 

PACS numbers: 75.10.Jm, 75.50.Ee 



Recent experimental studies of a spin-1/2 Kagome 
system ZnCu 3 (OH) 6 Cl 2 P, S, i show that the system 
is in a non-magnetic ground state. The Kagome lat- 
tice can be viewed as corner-sharing triangles in two- 
dimension(Fig[TJa)). The compound shows no magnetic 
order down to very low temperature (50mK) compared 
with the Curie- Weiss temperature (>200K). The spin 
susceptibility rises with decreasing temperature, but sat- 
urates to a finite value below 0.3 K. The specific heat is 
consistent with a linear T behavior below 0.5 K. There 
is no sign of a spin gap in dynamical neutron scattering. 
These observations led us to re-examine the issue of the 
ground state of the spin-1/2 Kagome lattice. 

Based on Monte Carlo studies of Gutzwiller projected 
wavefunctions, we propose the ground state to be a 
U(l)-Dirac spin-liquid state which has relativistic Dirac 
spinons. The low energy effective theory is a U(l) gauge 
field coupled to four flavors of two-component Dirac 
fermions in 2+1 dimension. This state was studied earlier 
in the mean- field approximation [4] . However, that study 
focused on an instability toward a Valence Bond Solid 
(VBS) state which breaks translation symmetry [4|; it was 
not appreciated that the U(l)-Dirac state can be a sta- 
ble phase. Using the Projective Symmetry Group 0,0,0] 
(PSG) technique, we reconsider the stability of the U(l)- 
Dirac state and find it can be stable. Our numerical cal- 
culations confirm that neighbor states like the VBS states 
and chiral spin-liquid state all have higher energies. 

One way to construct spin-liquid states is to introduce 
fermionic spinon operators d, Q /j and f± to represent 

the bosonic spin operator: Sj = 2Jia (Ja pfip- This re P" 
resentation enlarges the Hilbert space, and a local con- 
straint is needed to go back to the physical Hilbert space: 



model 



flfi = I- For the nearest neighbor Heisenberg 



H 



(i) 



can be decomposed via a Hubbard-Stratonovich trans- 
formation by introducing the complex field Xij living on 
the links. The path integral of the spin model is then 
Z = j dxdXdf 'dp ' e~ s , where the action is 



(2) 



we can substitute the spin operator by the spinon op- 
erators, so that the spin interaction is represented as a 
four-fermion interaction. The four-fermion interaction 



Here A is the Lagrangian multiplier to ensure the local 
constraint, and it can be viewed as the time component 
of a compact U(l) gauge field, whereas the phase of Xij 
can be viewed as the space components of the same gauge 
field. Only when the full gauge field fluctuations are in- 
cluded can one go back to the physical Hilbert space. 

With this fermionic representation, one can do a mean- 
field study of the spin- liquid states by taking Xij' s as 
mean-field parameters. For the Kagome lattice, the 
mean-field states are characterized by the fluxes through 
the triangles and the hexagons. Controlled mean-field 
studies were done by generalizing the SU(2) spin model 
to SU (N) spin model via introducing N/2 flavors of 
fermions [J, |lOj, and several candidate states were found: 

(i) Valence Bond Solid (VBS) states which breaks 
translation symmetry. 

(ii) a spin liquid state (SL-[^,0]) with a flux +7r/2 
through each triangle on Kagome lattice and zero- 
flux through the hexagons. This is a chiral spin 
liquid which breaks time-reversal symmetry. 

(iii) a spin liquid state (SL- [± ^ , 0] ) with staggered tt/2- 
flux through the triangles (+? through up trian- 
gles and —\ through down triangles) and zero- flux 
through the hexagons. 

(iv) a spin liquid state (SL-[5,7rl) with +7r/2-flux 
through the triangles and 7r-flux through the 
hexagons. 
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(v) a uniform RVB spin liquid state (SL-[0, 0]) with 
zero-flux through both triangles and hexagons. 
This state has a spinon Fermi surface. 

(vi) a U(l)-Dirac spin liquid state (SL-[0, 7r]) with zero- 
flux through the triangles and 7r-flux through the 
hexagons. This state has four flavors of two- 
component Dirac fermions. 

Marston and Zcng found that among the spin liquid 
states (jii|-(|vj), the chiral spin liquid SL- [f , 0] has the low- 
erest mean- field energy. But numerical calculations 11 
do not support a large chirality-chirality correlation, and 
Hastings |4[ found SL-[0, ir] to be the state with the low- 
est mean-field energy among the non-chiral spin liquid 
states. However its mean-field energy is still higher than 
that of All the above mean-field arguments are 

based on the expansion treatment of gauge fluctua- 
tion, which may fail when N — 2 in the physical case. 
To clarify which candidate is the physically low energy 
spin liquid state, we do a Monte Carlo study on the trial 
projected wavefunctions [l2j| . 

As we mentioned, fermionic representation enlarges 
the Hilbert space. One way to treat the unphysical 
states is to include gauge fluctuations. Another direct 
way to remove the unphysical states is to do a pro- 
jection by hand. Suppose that we have a mean-field 
ground state wavefunction ' m ean{x.ij)) with mean-field 
parameters Xij' s : we can project out all the unphysical 
states, then the resulting wavefunction \^ pr j) would be 
a physical state: \^ P rj(Xij)) = -Pol* rue an (Xij)), where 
Pd = — n i^ n i[) is the projection operator ensuring 
one fermion per site. It turns out that the calculation of 
energy (^ pr j\H\^ pr j) can be implemented by a Monte 
Carlo approach with power law complexity, which means 
that one can do a fairly large lattice. [l2| 

We note that states related by a global transformation 
Xij — * ~Xij represent the same spin wavefunction after 
projection. This is because projection of holes = f\ 
gives the same state as projection of fermions and the 
hopping parameter is transformed accordingly. This is a 
special case of the SU{2) gauge symmetry [l3j. In partic- 
ular, if Xij is real, the projected state and its energy are 
independent of a global sign change. 

For spin-1/2 nearest neighbor anti-ferromagnetic 
Heisenberg model Eq.|T|), we did the Monte Carlo cal- 
culation for energies of projected spin liquid states on 
lattices with 8x8 and 12x12 unit cells (each unit cell has 
3 sites). We choose mixed boundary conditions; i.e., pe- 
riodic along one Bravais lattice vector, and anti-periodic 
along the other Bravais lattice vector. The results are 
summarized in Table [T] 

We found that the U(l)-Dirac state (this state is the 
projected spin liquid of the mean-field state (jvjj proposed 
by Hastings[J|) has the lowest energy, which is — 0.429J 
per site. Note that these results change the order of 



Spin liquid 


8x8x3 lattice 


12x12x3 lattice 


SMS, ol 


-0.4010(1) 


-0.4010(1) 


sM±f,ol 


-0.3907(1) 


-0.3910(1) 


SL-[f,7r] 


-0.3814(1) 


-0.3822(1) 


SL-[0,0] 


-0.4115(1) 


-0.4121(1) 


SL-[0,tt] 


-0.42866(2) 


-0.42863(2) 



TABLE I: For all candidate projected spin-liquids, we list the 
energy per site in unit of J. The U(l)-Dirac state SL-[0, 7r] is 
the lowest energy state, and its energy is even lower than some 
numerical estimates of the ground state energy (see Table HI)) . 



Method 


energy per site 


Exact Diagonalizationfll] 


-0.43 


Coupled Cluster Method [14] 


-0.4252 


Spin-wave Variational method[15] 


-0.419 



TABLE II: We list the previous estimates for ground state 
energy in unit of J. 



mean- field energies of the spin liquids (fiT|) - (TvT|) . where the 
chiral spin liquid (JTTJ) was found to be of the lowest en- 
ergy. In Table [Til we list the estimates of the ground 
state energy by various methods. It is striking that even 
though the projected U(l)-Dirac state has no variational 
parameter, it has an energy which is even lower than some 
numerical estimates of ground state energy. Furthermore 
its energy is very close to the exact diagonalization result 
when extrapolated to large sample size. Thus we propose 
it to be the ground state of the spin-1/2 nearest neighbor 
Heisenberg model on the Kagome lattice. 

Hastings^ proposed a neighboring VBS ordered state 
as an instability of the U(l)-Dirac state. This state can 
be obtained by giving the fermions non-chiral masses. In 
particular, he proposed a VBS state with with a 2 x 2 
expansion of the unit cell. The 12 hopping parameters on 
the boundary of the star of David (six triangles suround- 
ing the hexagon) have amplitude xi> while all other hop- 
pings have amplitude \i- Our numerical calculations 
show that this VBS ordered state has higher energy (see 
Table ITlTI) . so the U(l)-Dirac state is stable against VBS 
ordering. Another neighbor state of the U(l)-Dirac spin 
liquid discussed by Hastings [H is obtained by giving the 
fermions chiral masses. The resulting state is a chiral 
spin liquid with broken time-reversal symmetry, and has 
6*-flux through the triangles and (tt — 2#)-flux through 
the hexagons (if 9 = the state goes back to the U(l)- 
Dirac state). In Table Hill we also show that non-zero 9 
increases the energy. 

To determine whether the U(l)-Dirac state is a stable 
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State 


8x8x3 lattice 


12x12x3 lattice 


U(l)-Dirac spin liquid 


-0.42866(2) 


-0.42863(2) 


VBS state(|xi/X2| = 1-05) 


-0.42848(2) 


-0.42844(2) 


VBS state(| X i/x2| = 0.95) 


-0.42846(2) 


-0.42846(2) 


Chiral spin liquid(0 = 0.05) 


-0.42857(2) 


-0.42853(2) 



TABLE III: We list the energy per site in unit of J for possible 
instabilities of the U(l)-Dirac spin liquid. The VBS state 
and the chiral spin liquid are those discussed in Ref[|( see 
text) . One can see that both VBS order and chiral spin liquid 
increase the energy. 



phase, we start with its effective theory 

+ ^ ~ vV'-ct] H , (3) 

IT 

where the first term comes from integrating out some 
higher energy fermions, and • • • represents other terms 
that are generated by interaction. The massless Dirac 
fermions in the effective theory come from the gap- 
less nodal spinons in the mean-field theory. The two- 
component Dirac spinor fields are denoted by ?p±cn where 
± label the two inequivalent nodes and a the up/down 
spins. Also, %p±„ = ?/4 ct t 3 , and the are Pauli ma- 
trices. The massless fermions lead to an algebraic spin 
liquid 0]. The stability of the U(l)-Dirac state can 
now be determined by examining the • • • terms: If • • • 
terms contain no relevant perturbations, then the U(l)- 
Dirac state can be stable. 

The potential relevant terms are 16 possible gauge- 
invariant, spin-singlet bilinears of ip± a . To see if those 
bilinears are generated by interaction or not, we need to 
study how lattice symmetries are realized in the effec- 
tive theory ((3]). Because spinons are not gauge invariant, 
lattice symmetry is realized in the effective theory as a 
projective symmetry, described by a PSG. This means 
that the realization of lattice symmetry includes nontriv- 
ial gauge transformations. We find that 15 of 16 bilinears 
are forbidden by translation symmetry and time-reversal 
alone. The remaining bilinear, which is allowed by sym- 
metry, is Yl± a i ) ±ai ) ±(T- This term shifts the spinon 
Fermi level to make the ground state to have exactly one 
spinon per site. In this case, the lower three of six spinon 
bands are filled and the spinon Fermi level is exactly at 
the gapless nodal points. This analysis tells us that the 
U(l)-Dirac state is stable in mean-field theory (and also 
in a large- N treatment). Because not all scaling expo- 
nents are known in such an algebraic spin liquid, pertur- 
bations other than fermion bilinears could in principle 
lead to an instability. However, so far, the variational 
wavefunction analysis suggests that this is not the case 
and that the U(l)-Dirac state is stable. 

Now we study the U(l)-Dirac spin liquid on the mean- 




(a) 

b 2 E(ky) 




(b) (c) 

FIG. 1: (a): We choose the six-site unit cell for the Kagome 
lattice. The sites are labeled 0, . . . , 5 as shown. Only those 
bonds depicted as bold solid lines and bold dashed lines are 
contained within the unit cell. The bold solid bonds have 
positive hopping, and the bold dashed bonds have negative 
hopping, (b): The Brillouin zone for the doubled unit cell 
(gray area), with reciprocal lattice basis vectors bi shown. 
The outer hexagon is the Brillouin zone for the 3-site unit cell 
of a single up-pointing triangle. The positions of the Dirac 
nodes are denoted by the black circles, (c): Plot of the band 
structure of the [/(l)-Dirac state on the line from fe = to 
fe = bi with energy in units of x>^(see text). The flat band is 
doubly-degenerate; all others are nondegenerate. The Fermi 
level corresponding to one spinon per site is indicated by the 
dashed line. 



field level. The U(l)-Dirac mean-field state is defined 
as the ground state of the following tight-binding spinon 

Hamiltonian: H mean = Jj2(ij) Xij fj a fi a + h - c - Alley's 
have the same magnitude and they produce zero-flux 
through the triangles and 7r-flux through the hexagons. 

Although the U(l)-Dirac state does not break transla- 
tion symmetry (because the translated state differs from 
the original state only by a gauge transformation), the 
unit cell has to be doubled to work out the mean-field 
spinon band structure. One can fix a gauge in which all 
hoppings are real as shown in Fig. [TJa). In this gauge 
the Dirac nodes are found to be at fc = (0, ±-^) as 
shown in Fig. QIb)(c), where a is the Kagome unit cell 
spacing, i.e., twice the nearest neighbor distance. These 
are isotropic Dirac nodes; i.e., the Fermi velocity is the 
same in all directions. More nodes can be obtained by 
covering the momentum space by repeating the Brillouin 
zones. One can easily see that the Dirac nodes form a 
triangular lattice in momentum space with lattice spac- 
ing -^2- . Note that the positions of the Dirac nodes are 
gauge dependent, but the momentum vectors connecting 
any two Dirac nodes are gauge invariant. Because the 
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spinon excitations are gapless at the nodal points, we ex- 
pect the spin-1 excitations of the U(l)-Dirac spin liquid 
are also gapless at zero momentum and those momenta 
connecting two Dirac nodes. 

For the ZnCu 3 (OH) 6 Cl 2 compound, the Heisenberg 
coupling was estimated to be J « 3OOK0, and one 
can calculate the Fermi velocity at mean-field level. We 
find that the Fermi velocity of the U(l)-Dirac ansatz is 
vf = where \ is the magnitude of the self-consistent 
mean-field parameter. Hastings^ found that x — 0.221. 
X describes the renormalization of the spinon bandwidth 
and is not expected to be given quantitatively by the 
mean-field theory. Hence in the formulae below we retain 
X as a parameter. We find vp = = 19% ■ 10 3 m/s. 

We can also calculate the specific heat of the U(l)- 
Dirac state at mean-field level. At low temperature 
(ksT <C xJ)i one expects a C oc T 2 law because of the 
Dirac nodes. The coefficient is related to vf- 



C _ 72C(3)7rfc|A 



T 2 



= l.l X ~ 2 ■ 1CT 3 Joule/mol K 3 , (4) 



where A is the area of the 2-D system. (Note that for 
ZnCu3(OH) 6 Cl2 compound, the unit cell spacing a = 
6.83A, so A = 2.4 ■ 10 5 m 2 /mol, where mole refers to one 
formula unit. We also used the fact that there are four 
two-component Dirac fermions.) 

In a magnetic field, the spinons will form a Fermi 
pocket whose radius is proportional to magnetic field 
strength. Therefore at low temperature ksT <C ^bB, 
the specific heat is linear in T: 



C _ 8tt kgAfisB 
T ~ 3(2Trhv F ) 2 



0.23 X ~ 2 B ■ 1CT 3 Joule/mol K 



where magnetic field B is in unit of Tesla. We also find 
in the temperature range \xbB <C ksT <C 



C 



24irAk 3 B T 2 - 



(2ithv F ) 2 



3C(3) 



21n2 /fi B B\ 2 
3 \k B TJ 



0(B 4 ) 



Keeping the lowest order correction, the specific heat has 
a temperature independent increase proportional to B 2 . 



16ir\n2k B A 2 
A ° = (2TTHVF) 2 iflBB) 

= 6.3x~ 2 B 2 ■ 10~ 5 Joule/mol K. 



(5) 



This is in contrast to the specific heat shift of a local mo- 
ment, which decreases with T as B 2 /T 2 . Eq.© provides 
a way to seperate the Dirac fermion contribution from 
that of impurities and phonons. 

The gauge field also gives a T 2 contribution to the 
specific heat. However, in a large- iV treatment this will 
be down by a factor of 1/N compared to the fermion 
contribution. Furthermore, the self energy correction due 
to gauge fluctuations does not lead to singular corrections 



to the Fermi velocity [l6j. so the T 2 dependence of C is a 
robust prediction. 

We notice that experiment observed that the specific 
heat of Kagome compound ZnCu3(OH)gCl2 behaves as 
C oc T 2 / 3 in zero magnetic field over the temperature 
window 106 mK < T < 600 rnKjU, which is enhanced 
from C oc T 2 law. This enhancement is suppressed by a 
modest magnetic fieldQ. Furthermore, over a large tem- 
perature range (10K to 100K), the spin susceptibility is 
consistent with Curie's law with 6% impurity local mo- 
ment. We propose that these impurity spins (possibly 
due to Cu located on the Zn sites) may be coupled to 
the spinons to form a Kondo type ground state with a 
Kondo temperature < IK, thus accounting for the large 
C/T and the saturation of the spin susceptibility below 
0.3K. The Kondo physics of impurities coupled to Dirac 
spinons is in itself a novel problem worthy of a seperate 
study. Meanwhile it appears to dominate the low tem- 
perature properties and obscure the true excitations of 
the Kagome system. We propose that a better place to 
look for the Dirac spectrum may be at higher temper- 
ature (above 10K) and as a function of magnetic field, 
where the impurity contributions may be suppressed and 
the unique signature of Eq. (j4} and Eq. ((5|) may be tested. 
On the other hand, we caution that from FigfjTc), the 
spinon spectrum deviates from linearity already at a rel- 
atively low energy scale(~ 0.5% J). Our theory also pre- 
dicts a linear T spin susceptibility of fc^T <C xJ ■ Knight 
shift measured by Cu NMR is the method of choice to 
seperate this from the impurity contribution. 

Finally we remark on a possible comparison with exact 
diagonalization studies which found a small spin gap of 
~ t^j and a large number of low energy singlets 11|. It 
is not clear whether these results can be reconciled with 
a [/(l)-Dirac spin liquid. Here we simply remark that 
in a finite system the Dirac nodes can easily produce a 
small triplet gap and that the gauge fluctuations may be 
responsible for low energy singlet excitations. 

We thank J. Helton and Y. S. Lee for helpful dis- 
cussions. This research is supported by NSF grant No. 
DMR-0433632 and DMR-0517222. 
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